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ABSTRACT 

In the common model supernova shock-acceleration of cosmic rays there are two open questions: 

1. where does the high energy cosmic rays below the knee (10^ - 10^ Gev) come from, and 

2. are cosmic ray accelerated only at their origin or contineuosly during their residence in 
the Galaxy. We show that 10^5 eV light nuclei are probably accelerted by associations of 



On ! supernovae. The ratio of the spectra of secondary to primary cosmic rays would be affected 
^ ' by repeated (or secondary) acceleration in the ISM during their propagation in the galaxy. 

The observed secondary and primary CR spectra are used to constrain the amount of such 
^ ' secondary acceleration by supernova remnants (SNR). Two cases are considered: weak shocks 



^ ■ (1 < M < 2) of old, dispersed remnants, and strong shocks (M > 3) of relatively young 
remnants. It is shown that weak shocks produce more secondary acceleration than what is 
permitted in the framework of the standard leaky box (SLB) model, making it inconsistent 
with dispersed acceleration that should be produced by SNR. If the SLB is modified to allow 
I a moderate amount of RA by week shocks, the RA produced by old SNRs agrees with the 
^ ' rate reqired to fit the secondary-to primaray cosmic-ray data, making a self consistent picture. 
Q\ • Significant secondary acceleration by strong shocks of young SNRs should lead to flattening of 
O 



the secondary-to primaray ratio at high energies, near ITeV/nucleon. 



5 ; WHAT IS SECONDARY ACCELERATION? 

'-^1 The cosmic ray spectrum in the range 1-10^ is probably produced by shock acceleration in the 
I ' Galaxy (e.g. Axford 1981). Its power-law shape J(p) ~ p^ '^ is beleived to be produced by strong 
2 ' shocks, probably of young SNR (Blandford and Ostriker, 1980), combined with escape from the 
c/3 . Galaxy. After the cosmic ray particle has left the SNR site and is propagating in the Galaxy, 
it may enconter another SNR, be trapped and be accelerated to a higher energy. However, the 
probability to encounter a SNR is proportional to its volume. It turns out to be significant 
^ I only for old SNR, hens with weak shocks. Indeed, Observations of the ratio of primary to 
d [ seconday cosmic rays suggest such a secondary acceleration process, but at the same time 
constrain its amount (Eichler 1980, Cowsick 1986). Wande et.al. (1987) have calculated the 
cosmic ray spectrum and the constraints on secondary acceleration set by the data under various 
secondary acceleration models, with the secondary acceleration amount and the shock-strength 
as parameters. Applying these models to RA by SNR (Wandel 1987; 1988) shows that one 
cannot avoid secondary acceleration if one assumes that the primary cosmic rays are produced 
by supernovae: the young SNR producing the the cosmic rays will expand and occupy a large 
enough fraction of the Galaxy to contribute a considerable amount of secondary acceleration . 
The standard leaky-box model has therefore to be modified to include RA by SNR. 

BEYOND THE KNEE 

Acceleration by supernova shocks can produce cosmic rays with energy up to 10^^ eV. It is 
difficult to reach higher energies, because when the gyroradius of the accelerated particles 
becomes of the order of the shock size, the shock-acceleration mechanism cannot function, as 
the particle escapes. The gyroradius is given by 



(1) 



so that the maximal energy is 

E^r, = 0.1R(pc)ZB(iiG). (2) 

where £'15 is the maximal energy in units of 10^^ eV. To reach energies beyond 10^^ eV, larger 
shocks are required, and those can be obtained by going to older SNR and to associations of 
supernovae. Older SNR have weaker shockes, which produce a steaper power law. Indeed, a 
significant steepening is observed in the cosmic ray spectrum beyond 10^^. This could indicate 
that old SNR associate to produce shocked volumes of sizes of lOOpc and more. With galactic 
magnetic fields of 3/iG lOOpc shocks, if abundant enough, could produce the observed cosmic 
ray spectrum up to 3 lO^^eV for protons and up to lO^^eV for iron. The formalism presented 
by Wandel (1988) is used to convolve the size of the SNR with teir age and shock strength, to 
determine the probability and rate of secondary acceleration , and the efffective slope of the 
produced cosmic ray spectrum. 

THE SECONDARY ACCELERATION MODEL 



The steady state distribution of a particle population subject to contineous acceleration 
during its propagation is described by the integral equation 
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where Jo is the source distribution,_R = Roip/Po) " and S are loss terms due to escape from 
the galaxy and spallation, B is the rate of encountering accelerating events and 
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is the index of the power law (J ~ p~'^) produced by the shock of Mach number M (q=2 for 
strong shocks, and larger for weaker shoks). The distribution of each secondary species satisfies 
a similar equation, with the source term replaced by the sum over all primary species Jj which 
contribute to the secondary in question, Js{p) = T^SigJiip), with Sis is the spallation rate from 

i to s, and Jj satisfies equation (3). 



SECONDARY ACCELERATION AND THE LEAKY BOX MODEL 



In the framework of the Standard Leaky Box (SLB) model no secondary acceleration is 
allowed, except the margin due to the observational uncertainties. If the SLB model is modified 
by allowing a modest amount of secondary acceleration (and assuming a larger escape rate, 
which compensates for the increase in the secondary to primary ratio caused by secondary 
acceleration ), the constraint on secondary acceleration is less stringent. Wandel et.al. (1987) 
have calculated the secondary to primary (boron/carbon, denoted by B/C) energy spectrum 
for a modified leaky box models (MLB) with various escape laws, secondary acceleration rates 
and shock indices q. Their best fit of the B/C data permits a significantly larger secondary 
acceleration rate than in the SLB model. For shocks with q — A (M = 1.7) 

B < 0.03 Ro = 0.11 (SLB) 

B ^ (0.2 ± 0.05) Ro = 0.2 (MLB) 

where B and i?o are measured in units of inverse path length (gcm"^) 

Similar constraints can be derived from the primary spectra of protons and alpha particles 
(e.g. Webber, Golden and Stephens 1987). 
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ACCELERATION BY SUPERNOVA REMNANTS 
Weak shocks 



The number of supernova remnants encountered by a cosmic ray particle during its res- 
idence time in the galactic disk, (r^ ~ lO'^yr at a few GeV/nucleon) is given by Ncr = 
STcr{^+CY , where S is the supernova rate (per unit volume) in the galaxy, and C = rcr{t)/r 
represents the cosmic ray propagation during the remnant's lifetime. The average cosmic 
ray residence time in the galactic disk can be written in the form T^r — {l.ScmpnRo)"^ ~ 
510^i?Q ^yr. The reacceleration rate due to supernova remnants becomes 

3 

^O.USi2n~%i{l + C,)\ (6) 

where u = dr/dt = 100u2km/s is the shock velocity, the shock index q is related to the 
Mach number M = u/cg by q{M) and 5*12 = 5'/10~^^ pc~^yr~^ (one supernova every 30 
years corresponds to 5'i2 = 0.5). We take the sound speed in the ISM = 150/cm/s, so 
M — l.l{q — 4) corresponds to it = 250km/s. If the effective density seen by supernova 
remnants is that of the warm component of the ISM, n^ff = 0.1 (McKee and Ostriker 1977), 
then ~ 23pc, Cq ~ 0.2 — 0.4 (depending on the cosmic ray propagation, Wandel 1988), and 
eq. (6) gives B = 0.2 — 0.4. While this figure is by a factor 10 larger than the constraint in the 
SLB model, it is just the secondary acceleration rate required by the modified model in order 
to fit the B/C data. 

Strong-Shock Secondary Acceleration by Young SNRs 



B{q) = N,rRo = 3.0 Surl^ 
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The effect of secondary acceleration by young SNRs on the secondary to primary ratio 
has been considered in previous papers (Wandel 1990; 1991). Although the probability of 
encountering a young SNR is small, its effect on the cosmic ray distribution at high energies 
is strong, and may be detected. At high energies the solutions for the cosmic ray distribution 
(eq. 3) become assymptotically power laws, with the power index depending on the source 
acceleration index qo, the secondary acceleration index q, and the escape law R oc p'"'. It can 
be shown that the secondary to primary ratio has the assymptotic form (Wandel et.al. 1987) 

^ r p-a q>qo + a . . 

Jp q<qo + a ^ ^ 

In particular, for q = qo the secondary to primary ratio should fiatten at high energies. The 
energy at which this fiattening occurs depends on the secondary acceleration rate B and on the 
secondary acceleration index q. Prom the observed cosmic ray spectra we have qo + a 2.7 
and q; ~ 0.5 — 0.6. A significant fiattening of the B/C ratio requires shocks with q < 2.4, which 
gives M > 4 and u = 600km/ s. Young supernova remnants may evaporate the cold phase of 
the ISM (McKee and Ostriker 1977), in which case rieff = 1. If this is not the case, only the 
warm component is swept by the shock, so rie/f = 0.1. Pinally, If the supernova expands within 
a cavity of hot ISM, Ueff = 0.003 - 0.01. Eq. (6) gives 



B{q < 2.4) 



0.003 n = 1 (evaporation) 
0.03 n = 0.1 (no evaporation) 
0.3 n = 0.01 (hot ISM) 



If the transition between evaporative and non-evaporative expansion happens at an intermediate 
velocity (Cowie, McKee and Ostriker 1981), B will be between 0.003 and 0.03. The secondary to 



primary (Boron/Carbon) ratio has been observed up to ~1 TeV by the Spacelab-2 experiment 
of the University of Chicago (Meyer et.al.l987). The poor statistics at the high energy end 
prevent a decisive conclusion. However, prehminary data suggest that the B/C ratio in the 
0.1-lTcV range is approximately 0.1. Comparing this with the theoretical calculation we find 
that such a B/C ratio requires a secondary acceleration rate B{q = 2 — 3) ~ 0.03 — 0.1, so 
strong-shock secondary acceleration is probably important, and young supernova remnants see 
a relatively low effective density. Even if the B/C ratio beyond lOOGeV/n is eventually found to 
decrease, strong shock secondary acceleration by evaporative young supernova remnants places 
a lower limit on the secondary acceleration rate at B — 0.003, which would predict a flattening 
near ITeV. 
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